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Abstract 



Processes of collimation of electrically charged particles near a rotating black hole 
are discussed. It is assumed that the black hole is immersed in a weak magnetic 
field aligned with rotation axis. This situation is relevant for understanding pre- 
collimation of astrophysical jets. Magnetic field affects motion of material and re- 
stricts validity of various scenarios which adopt the test-particle (cold plasma) ap- 
proximation. A simplified criterion which estimates relevance of such approximation 
is discussed in connection with the mechanism of the dissipative collimation, as pro- 
posed by de Felice & Curir [19]. 
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1 INTRODUCTION 

Numerous observations confirm that collimated outflows of matter are rather 
generic phenomenon connected with certain types of astronomical objects [27]. 
These jets exist on different length-scales and they are associated with various 
types of sources ranging from stars to galactic nuclei — i.e. over nine orders of 
magnitude in the mass of the central source. The origin of jets is probably di- 
verse but they share common properties. For example, it has been speculated 
about analogies between electromagnetic processes which accelerate particles 
near pulsars [5,39] versus processes in magnetospheres of supermassive black 
holes in active galactic nuclei (AGN) [10,42]. AGNs with jets show many di- 
verse properties, but it is tempting to link the differences, at least partially, 
to orientation of these objects with respect to the observer. The interest in 
this subject has been amplified by recent discoveries of relativistic outflows in 
our Galaxy [26,40,59] which have their well-known counterpart in extragalac- 
tic superluminal jets [65]. Unifying schemes have been proposed for cosmic 
sources with jets. Useful review articles summarizing our knowledge can be 
found in the literature, both for stellar-scale objects in the Galaxy [45] and 
for extragalactic jets [4,60]. 

Extragalactic jets are presumably formed in the innermost regions of the 
source (within a few or a few tens of gravitational radii, R g , from the center) 
and they emanate outwards along rotation axis of the central object. We as- 
sume the model with a compact rotating object in the center. Gravitational 
field of the central object is approximated by the Kerr metric while self-gravity 
of the jet material and accretion disk is neglected. It has been widely recog- 
nized that large-scale outflows can be adequately described in the magnetohy- 
drodynamic regime but situation very close to the horizon is less understood. 
The initial phase of the jet formation is sometimes called pre- collimation, in 
distinction to the processes of successive collimation which operate in more 
distant regions. 

Several mechanisms of the jet pre-collimation have been proposed. Mag- 
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netic fields play most probably a major role in focusing and maintaining col- 
limated outflows on their course. Numerous authors have studied collimation 
within the hydromagnetic framework [9,12,13,37,44]. It has been shown that 
toroidal component of the magnetic field is maintained by rotation of the 
accretion disk and enhances collimation [61]. We were interested in the contri- 
bution of the Kerr geometry to the resulting collimation of jets, and it appears 
quite natural to presume that our attention can be restricted to the region 
within a few tens of gravitational radii from the center. The inner region of 
the object is crucial for the theory of jet formation but it remains beyond 
current observational capabilities which are of the order of (10 2 -10 3 ) R g for 
extragalactic sources (the best linear resolution, approximately 0.01 pc, has 
been achieved with the radio jet in the galaxy 3C274; [28]). 

Abramowicz & Piran [1] and Sikora & Wilson [57] consider collimation 
inside a funnel of a luminous thick accretion disk. Material of the jet is in mu- 
tual interaction with the disk radiation which determines its terminal speed 
[48]. Nowadays, the idea of extremely thick disks with very narrow funnels 
and highly super-Eddington luminosity (e.g., [21]) is not favoured because 
these models suffer from several inconsistencies but the general scheme of jets 
flowing along the disk axis remains viable with more sophisticated models of 
accretion disks. The model has been advanced by detailed quantitative calcu- 
lations of acceleration/deceleration of jets due to radiation pressure and losses 
due to cooling, both within the framework of the hydrodynamic [16,36,47] 
and the test-particle approximation [38,56,57,62]. These calculations impose 
strong limitations on radiatively driven jets because their material cannot 
reach Lorentz factors much greater than unity. Hydromagnetic scenario is 
thus currently favoured [3,44]. 

The above-mentioned schemes consider a massive compact object to be 
present in the center of the source but the Kerr geometry and the effects of 
general relativity are not crucial for their function. On the other hand, de Fe- 
lice & Curir [19] and de Felice & Carlotto [18] have explored a special family of 
geodesies spiralling along the axis of the Kerr black hole (vortical geodesies; see 
[17,20]) and determined constraints on the rate of change of energy and angu- 
lar momentum of outflowing material that may result in collimation. They did 
find collimation but the physical nature of dissipative processes that cause the 
loss of particle energy and angular momentum remains unclear. It is the aim 
of the present contribution to advance the latter model by systematic discus- 
sion of the parameter space of particle trajectories. We want to determine, by 
a simplified but systematic approach, circumstances under which the above- 
mentioned scheme could be relevant (rather than build up our own model). 
In other words, we ask whether specific features of motion in Kerr geometry 
are relevant for pre-collimation of astrophysical jets. Indeed, the original mo- 
tivation to deal with this problem was a suspicion that even relatively weak 
magnetic fields impose strong limitations upon the model. (This fact has been 
quoted and applied in numerous works; we wish to discuss the problem in a 
more systematic way with the Kerr geometry.) Our arguments could be ap- 
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plied also to other models based on the properties of the geodesic motion [7], 
spin-curvature coupling (another purely general relativistic effect, [54]), and in- 
fluence of magnetic fields which affect the spacetime geometry [30]. The latter 
process requires dimensionless product (3 of magnetic intensity and the mass of 
the central object to be of the order of unity, which is value much higher than 
= 5 x 10- 8 B 4 M 8 « fCT 8 , where B 4 = 5/(f0 4 gauss) and M 8 = M/(f0 8 M ) 
denote typical values relevant for the innermost magnetosphere in AGNs. Di- 
rect interaction of the jet material with external electromagnetic fields appears 
more important. (Geometric units will be used hereafter, c — G — 1; time, 
mass, electric charge and [magnetic intensity] _1 have dimension of length, i.e. 
cm. Next, quantities with dimension of length will be expressed in units of 
M.) 

Structure of this text is as follows: First, equations for generalized energy 
and generalized angular momentum of electrically charged particles are given. 
Dimensionless parameters characterizing, locally, electromagnetic effects upon 
the test-particle motion are then defined and evaluated. It should be noted 
that in the present contribution we do not address the problem of extraction of 
rotational energy from a magnetized black hole which also relies on properties 
of the Kerr spacetime [25,43,58,64], neither we discuss purely electromagnetic 
collimation due to toroidal magnetic fields. 

2 CHARACTERISTIC LENGTH-SCALE 

2. 1 Details of the model 

We assumed that a weak magnetic field is generated far from the black hole, 
outside the region where motion of matter is studied. The spacetime is de- 
scribed by the Kerr metric with an asymptotically uniform magnetic test field 
[63] . This solution reflects the large-scale field which is generated far from the 
black hole. We further assumed that the field is aligned with rotation axis 
because non-aligned fields exert torque on the black hole, so that its rotation 
slows down and the black hole gets aligned with the field [33,49,51]. Typical 
time-scale for the black-hole alignment due to electromagnetic torques is pro- 
portional to (3~ 2 and rather long (in physical units, r/(10 10 yr) ps M% l B~[ 2 ) 
but tidal effects enhance alignment and decrease r [34,52]. 

Structure of an asymptotically uniform magnetic field is simple and colli- 
mation by Kerr geometry is presumably easier to notice [53]. (As mentioned 
above, toroidal fields, B^, contribute to purely electromagnetic collimation 
and to the total power-output from a black hole [46]; in time-dependent cal- 
culations, B^ = is often taken as an initial configuration of the magnetic 
field [44]). Jets are formed near a black hole but the material which forms jets 
is an open question (see, e.g., [3,35] for general issues related to formation of 
jets). In the case of electron-positron plasma, the specific electric charge of in- 
dividual particles is \q\ ss 2 x 10 21 (elementary charge q and mass of electron m 
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will be assumed in numerical estimates, q = q/m). Dimensionless parameter 
e = \/3q\ determines motion of charged particles. 

Trajectory of individual particles is determined by Larmor gyrations around 
magnetic lines of force (Larmor radius r L ) and the drift motion across the 
field-lines (e.g., [2]). Motion of electrically charged particles near a black hole 
has been studied by numerous authors [8,15,22,55]. Prasanna [50] and Karas 
& Vokrouhlicky [30] show examples of test-particle trajectories near magne- 
tized black holes. Now we are interested in the bulk motion of the material 
rather than individual trajectories which can be very complicated. It is thus 
useful to disregard local gyrations and accept the guiding-center approxima- 
tion. Damour et al. [14] determined shape of the plasma flow-lines near a 
weakly magnetized Kerr black hole in the guiding-center approximation. Dis- 
cussion has been generalized to the case of an electrically charged rotating 
black hole by Hanni & Valdarnini [24] (with a weak magnetic field), and Karas 
& Vokrouhlicky [31] (an exact solution of the Einstein-Maxwell equations with 
an arbitrarily strong magnetic field). These authors verified that plasma moves 
along rotation axis and they concluded that asymptotically uniform magnetic 
field evidently contributes to collimation. 

The situation becomes much less understood when mechanism proposed by 
de Felice & Curir [19] (motion along vortical geodesies plus tiny, yet undeter- 
mined dissipation of energy and angular momentum) is taken into account. 
One can, however, estimate relevance of this scheme by calculating the relative 
change of orbital parameters due to magnetic field. 

2.2 Rate of change of orbital parameters 

In this section, we will estimate the rate of change of orbital parameters (en- 
ergy and angular momentum with respect to rotation axis) of individual par- 
ticles and integrate the rate over particle distribution. This approach gives us 
a local criterion of importance of electromagnetic forces. 

The model is described by an axisymmetric stationary spacetime metric g^ u 
(/i, z/ = 0, . . . , 3). Electromagnetic test field is characterized by the tensor of 
electromagnetic field, and by corresponding four-potential: = A v ^ — A^ u . 
It should be emphasized that we do not have to impose further restrictions 
upon gravitational and electromagnetic fields, apart from the above mentioned 
assumption about axial symmetry and stationarity. To be specific, however, 
we will consider the Kerr metric [41] with an aligned asymptotically uniform 
magnetic field. Structure of the magnetic field has been explored by several 
authors (e.g., [23,49] in the case of aligned fields, and [6,32,29] in the case of 
inclined fields). In Boyer-Lindquist coordinates, x M = {t,r, 9, 0}, 





= (3 sin 2 9 [\ (r 2 + a 2 ) - aVlT 1 (l + cos 2 0) . 



(2) 
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Here, E = r 2 + a 2 cos 2 9. Each timelike geodesic (four-momentum p* 1 = m-u' 1 ) 
is associated with two conserved quantities — specific energy 



^0 = —U t \ q =o, (3) 

and specific angular momentum with respect to the symmetry axis 

A =M^|g = . (4) 

Four-momentum along a trajectory of an electrically charged particle is deter- 
mined by equation 

and corresponding conserved quantities are: 

e = - (u t + qA t ) , X = u 4> + qA^ (6) 

(generalized energy and angular momentum component, respectively). 

Any scenario which explains collimation of particles in terms of curvature 
effects acting upon free particles (as in Bicak, Semerak & Hadrava [7]) or par- 
ticles that are slightly perturbed by dissipative forces (as in de Felice et al. 
[18,19]) turns out to be astrophysically irrelevant when electromagnetic forces 
strongly affect motion of particles. In order to introduce a quantitative cri- 
terion for electromagnetic effects, one can define two parameters, Se(r, #;p M ) 
and SX(r,6;p^), which characterize the relative change of eo and Ao: 

, £o\xv+ u vdr _ . F tr u r + F te u e . . 

1 ■ = q A — da = oe da, (7) 

e 1 a;^ qA t + e 

1 - A °' wd - = qK F **" r+ _ F ;+ U ° da = 5\ da. (8) 
Ao | if A — qAff, 

Here, da denotes interval of proper time scaled with the light-crossing time 
across the characteristic length-scale A: da = dr/A. 

It is further postulated that A oc r L = Tvq' 1 ^ 1 (with V = 1/^l-v 2 < 10 
denoting the Lorentz factor, v velocity with respect to the locally non-rotating 
frame). Hence, (Se) and (SX) become independent of magnetic field strength. 
This fact is understandable: the Larmor radius decreases with (3 increasing, 
which means that the characteristic length along which the change of pa- 
rameters is determined in eqs. (7)-(8) decreases as well. Such a choice of A 
is well-founded (for our purpose of an order-of-magnitude criterion) because 
geodesic trajectory certainly cannot approximate a real trajectory of a charged 
particle on a scale greater then r^. (On the other hand, it is easy to evaluate 
(Se), (SX) also for a different choice of A; (3 then becomes another free param- 
eter, however.) Considering the collimation processes acting on length-scales 
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Fig. 1. Graphs of (Se) [panels (a), (c)] and (SX) [panels (b), (d)] as a function of 
radius x = 1 - R g /r (R g = 1 + \/I - a 2 , < x < 1) and (0 < 9 < tt). Upper 
panels, (a)-(b), show mean values taken over all trajectories with T < r max while 
lower panels, (c)-(d), deal with vortical trajectories only (notice a narrow gap near 
= 7r/2: vortical trajectories do not cross equatorial plane). Here, a = 1, < Y < 3, 
s = 0. See the text for details. 



Fig. 2. As in Fig. 1 but for < T < 2, s = 2. 

of a few (as in de Felice & Curir [19]), A should be comparable with or 
greater than R g ; this is further restriction on the upper limit for (3. 

Both parameters, Se and 5\, are denned locally (i.e. r, 9 given) and they also 
depend on particle's p M . We study the region outside the black-hole horizon 
by averaging over distribution of particles in the momentum space. We define 

(Se) 2 = —jdn J drn(r)fe 2 , (9) 

4?r r min 
and analogously 

r ma x 

(S\) 2 = l-JdQ J drn(r)5A 2 . (10) 

4^ r min 

Integration is taken over the particle distribution in energy, n(T), and over 
all directions of their local velocity. Values of (Se) > 1, (SX) > 1 mean that 
approximation of geodesic motion with a small perturbation is inappropriate, 
while values of both parameters much less than unity, simultaneously with 
A > R g , indicate that this approximation might be meaningful. 



3 RESULTS 

We evaluated parameters (Se) and (SX) for the Kerr metric and electromag- 
netic test field (l)-(2). Equations (3)-(4) can be written in the explicit form: 

e = VSAA" 1 T + LUX , (11) 

A = VAS- 1 Tv 4, sin 2 9, (12) 

where A = r 2 — 2r + a 2 , A = (r 2 + a 2 ) 2 — Aa 2 sin 2 9, u = 2arA~ x are functions 
from the Kerr metric in standard notation [41]; is azimuthal component of 
the speed of particle with respect to LNRF. Energy distribution of particles 
was approximated by a power-law: n(T) oc T~ s in a restricted interval of 
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energy (T max < 20, < s ^ 2; power-law energy distribution is motivated by 
astrophysically relevant situations). 

Figures 1-2 illustrate our results for an isotropic distribution of particles 
with respect to LNRF. We assumed A xs r L (characteristic length) for defi- 
niteness. In these two figures, (Se) <C 1, (SX) <C 1 [panels (a) and (b)], and 
one concludes that approximation of the close-to-geodesic motion may be rele- 
vant for modelling the jet precollimation in the given region of (r, 9). We have 
also examined separately the case of vortical trajectories which play a crucial 
role in discussion of de Felice & Curir [19] [panels (c) and (d)]. In the latter 
case, the assumption about isotropic distribution is supplemented by specific 
conditions for vortical trajectories: 

f>0, |L|<a 2 f, L<\ 2 <^±£^ (13) 

4a z T 

(L denotes the fourth constant of the Kerr metric, T = e\ — 1). Regions in 
the (r, 6) plane with (8e) ^ 1, (SX) ^ 1 have been excluded from Figures for 
clarity (particularly, close to the horizon and rotation axis). 

Figures 1-2 represent a typical situation. Further illustrations in which the 
parameter space is investigated systematically can be found on World-Wide 
Web. 1 



4 CONCLUSION 

We derived an order-of-magnitude criterion for possible relevance of those 
models of the jet pre-collimation (at distances of few R g from the black hole) 
which are based upon test-particle approximation: (Se) 2 , (SX) 2 <C 1 (inde- 
pendently of B), and A > R g (depends on the value of B: A oc B^ 1 ). Our 
estimates employ characteristic length-scales which can be shorter when col- 
lisions of particles in plasma are important or if the magnetic field is domi- 
nated by a short-scale chaotic component. This means the region of magnetic 
dominance can be larger than our criterion indicates. Although only rough 
estimates of strength of the magnetic field are currently available, it appears 
that the guiding-center approximation is the most relevant approach among 
the models based on test-particle motion. One should note that general rela- 
tivity effects remain important for motion of material not only because of the 
presence of the black hole in the center but also because the structure of elec- 
tromagnetic field itself is affected by strong gravity. It is also worth noticing 
that explicit formulae for non-aligned test fields are known and can be studied 
in a similar way. 

Our present discussion has been restricted by assumptions about the large- 
scale structure of the magnetic field, isotropic distribution of particle velocities 
in LNRF, and power-law distribution of their energy. We do not expect our 

1 http:/ /otokar. troja.mff.cuni.cz/user/karas/au_www/karas/papers. htm 
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results to be sensitive to these assumptions unless the particle distribution is 
very anisotropic; it should be repeated that farther from the center hydromag- 
netic description is appropriate. 

ACKNOWLEDGEMENT 

The authors acknowledge helpful comments and critical remarks from pro- 
fessor Fernando de Felice, and from participants of the Seminar on General 
Relativity in Prague. We acknowledge suggestions by an anonymous referee 
concerning definition of A. We thank Joan Evans for reading the manuscript. 
This work has been supported by the grants GACR 205/97/1165 and GACR 
202/96/0206 in the Czech Republic. 

References 

[1] Abramowicz, M. A., Piran, T. (1980). Astrophys. J., 241, L7. 
[2] Alfven, H., Falthammar, C.-G. (1963). Cosmical Electrodynamics (Clarendon 
Press, Oxford). 

[3] Begelman, M. C. (1995). Proc. Natl. Acad. Sci., 92, 11442. 
[4] Begelman, M. C, Blandford, R. D., Rees, M. J. (1984). Rev. Mod. Phys., 
56, 255. 

[5] Benford, G. (1984). Astrophys. J., 282, 154. 

[6] Bicak, J., Janis, V. (1985). Mon. Not. R. Astr. Soc, 212, 899. 

[7] Bicak, J., Semerak, O., Hadrava, P. (1993). Mon. Not. R. Astr. Soc, 263, 545. 

[8] Bicak, J., Stuchh'k, Z., Balek, V. (1989). Bull. Astron. Inst. Czechosl, 40, 65. 

[9] Blandford, R. D., Payne, D. G. (1982). Mon. Not. R. Astr. Soc, 199, 883. 
[10] Blandford, R. D., Znajek, R. L. (1977). Mon. Not. R. Astr. Soc, 179, 433. 
[11] Bridle, A. H. (1985). Can. J. Phys., 64, 353. 
[12] Camenzind, M. (1986), Astron. Astrophys., 162, 32. 

[13] Chakrabarti, S. K., Bhaskaran, P. (1992). Mon. Not. R. Astr. Soc, 255, 255. 
[14] Damour, T., Hanni, R. S., Ruffini, R., Wilson, J. R. (1978). Phys. Rev. D, 
17, 1518. 

[15] Dymnikova, I. G. (1986). Usp. Fix. Nauk, 148, 393. 

[16] Falcke, H., Biermann, P. L. (1995). Astron. Astrophys., 293, 665. 

[17] de Felice, F., Calvani, M. (1972). II Nuovo Cimento, 10B, 447. 

[18] de Felice, F., Carlotto, L. (1995). Preprint PD 95/th/49. 

[19] de Felice, F., Curir, A. (1992). Class. Quant. Grav., 9, 1303. 

[20] de Felice, F., Clarke, C. J. S. (1990). Relativity on curved manifolds (Cambridge 

University Press, Cambridge). 
[21] Frank, J., King, A. R., Raine, D. J. (1995). Accretion Power in Astrophysics 

(Cambridge University Press, Cambridge). 
[22] Gal'tsov, D. V., Petukhov, V. I. (1978). Zh. Eksp. i Teor. Fiz., 74, 801. 
[23] Hanni, R. S., Ruffini, R. (1976). Lett. Nuovo Cimento, 15, 189. 
[24] Hanni, R. S., Valdarnini, R. (1979). Phys. Lett., 70A, 92. 
[25] Hirotani, K., Takahashi, M., Nitta, S.-Y., Tomimatsu, A. (1992). Astrophys. J., 

386, 455. 



9 



Hjellming, R. M., Rupen, M. P. (1995). Nature, 375, 464. 

Hughes, P. A. (ed) (1991). Beams and Jets in Astrophysics (Cambridge 

University Press, Cambridge). 

Junor, W., Biretta, J. A. (1995). Astron. J., 109, 500. 
Karas, V. (1989). Phys. Rev. D, 40, 2121. 

Karas, V., Vokrouhlicky, D. (1990). Class. Quant. Grav., 7, 391. 
Karas, V., Vokrouhlicky, D. (1991). J. Phys. I (France), 1, 1005. 
King, A. R. (1977). Math. Proc. Camb. Phil. Soc, 81, 149. 
King, A. R., Lasota, J. P. (1977). Astron. Astrophys., 58, 175. 
Kumar, S., Pringle, J. E. (1985). Mon. Not. R. Astr. Soc, 213, 435. 
Konigl, A. (1985). Can. J. Phys., 64, 362. 

Li Zhi-Yun, Chiueh Tzihong, Begelman, M. C. (1992). Astrophys. J., 394, 459. 
Lovelace, R. V. E., Mehanian, C, Mobarry, C. M., Sulkanen, M. E. (1986). 
Astrophys. J.S, 62, 1. 

Melia, F., K6nigl, A. (1989). Astrophys. J., 340, 162. 

Michel, F. C. (1982). Rev. Mod. Phys., 54, 1. 

Mirabel, I. F., Rodriguez, L. F. (1994). Nature, 371, 46. 

Misner, C. W., Thorne, K. S., Wheeler, J. A. (1973). Gravitation (W. H. 
Freeman and Co., New York). 

Nitta, S., Takahashi, M., Tomimatsu, A. (1991). Phys. Rev. D, U, 2295. 

Okamoto, I. (1992). Mon. Not. R. Astr. Soc, 254, 192. 

Ouyed, R., Pudritz, R. E., Stone, J. M. (1997). Nature, 385, 409. 

Padman, R., Lasenby, A. N., Green, D. A. (1991). In Hughes, P. A. (ed) Beams 

and Jets in Astrophysics (Cambridge University Press, Cambridge) p. 484. 

Park, S. J., Vishniac, E. T. (1994). Astrophys. J., 426, 131. 

Peter, W., Eichler, D. (1996). Astrophys. J., 466, 840. 

Phinney, E. S. (1987). In Zensus J. A., Pearson, T. J. (eds) Superluminal Radio 

Sources (Cambridge University Press, Cambridge) p. 301. 

Pollock, M. D., Brinkmann W. P. (1977). Proc. R. Soc. Bond. A, 356, 351. 

Prasanna, A. S. (1980). Rivista del Nuovo Cimento, 3, 1. 

Press, W. H. (1972). Astrophys. J., 175, 243. 

Scheuer, P. A. C, Feiler, R. (1996). Astrophys. J., 282, 291. 

Semerak, O. (1995). Czech. J. Phys, 45, 1. 

Semerak, O. (1997). Private communication (to be submitted). 

Sharp, N. A. (1979). Gen. Rel. Grav., 10, 659. 

Sikora, M., Sol, H., Begelman, M. C, Madejski, G. M. (1996). Mon. Not. R. 
Astr. Soc, 280, 781. 

Sikora, M., Wilson, D. B. (1981). Mon. Not. R. Astr. Soc, 197, 529. 
Takahashi, M., Nitta, S., Tatematsu, Y., Tomimatsu, A. (1990). Astrophys. J., 
363, 206. 

Tingay, S. J., Jauncey, D. L., Preston, R. A., Reynolds, J. E., Meier, D. L. et al. 
(1995). Nature, 374, 141. 

Urry, C. M., Padovani, P. (1995). Publications of the Astronomical Society of 
Japan, 107, 803. 



10 



[61] Ustyugova, G. V., Koldoba, A. V., Romanova, M. M., Chechetkin, V. M., 

Lovelace, R. V. E. (1995). Astrophys. J., 439, L39. 
[62] Vokrouhlicky, D., Karas, V. (1991). Astron. Astrophys., 252, 835. 
[63] Wald, R. M. (1974). Phys. Rev. D, 10, 1680. 
[64] Wagh, S. M., Dadhich, N. (1989). Phys. Rep., 183, 137. 

[65] Zensus, J. A., Pearson, T. J. (eds) (1990). Superluminal Radio Sources 
(Cambridge University Press, Cambridge). 



11 



